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ABSTRACT 


The Explorer 35 spacecraft experiences a small, temporary secu- 
lar decrease in spin period (approximately 2 parts in 1000) during every 
pass through the lunar shadow. Because of this, the pseudo- sun pulse, 
generated when the spacecraft is in the shadow, no longer accurately 
indicates the sunward direction. This report describes an empirical 
method of fitting a spin period of the form a + b exp(-t/t 0 ) to the as- 
pect information available before and after the optical shadow (t is 
time since the spacecraft entered the shadow, and a, b, and t 0 are 
constants to be determined). The resulting aspect information is esti- 
mated to be within ±10 degrees of the true direction. 
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ASPECT DETERMINATION IN LUNAR 
SHADOW ON EXPLORER 35 


by 

Harold E. Taylor 

Goddard Space Flight Center 


INTRODUCTION 

Explorer 35, injected into selenocentric orbit on July 22, 1967, experiences a decrease in its 
moment of inertia as it passes through the lunar shadow. This apparently is caused by the cooling 
of the spacecraft and, most probably, by inward contraction of the solar paddles. Because angular 
momentum is conserved, the spin period decreases slightly (a few parts per thousand), and the 
pseudo-sun pulse generated onboard the spacecraft no longer accurately indicates the direction of 
the sun. This report describes an empirical approach to solving the resulting satellite aspect 
problem. 


STATEMENT OF THE PROBLEM 

A variable spin period of the form 

t = a + be " t / t o (1) 

is assumed, during the passage through shadow 
(a, b, and t 0 being fitted parameters to be de- 
termined). The spin periods (which are meas- 
ured during the recovery following a passage 
through shadow) obey this type of exponential 
law (Figure 1). During the shadow, the absence 
of a true sun pulse precludes the accurate de- 
termination of the spin period. However, to the 
extent that the azimuth angle of the magnetic 
field vector in the satellite spin plane remains 
constant throughout the shadow, the magnetic 
field measurements themselves give an esti- 
mate of the satellite azimuth. A comparison of 



Figure 1 -Observed r - a versus time for the recovery 
of the spacecraft following the first shadow pass for 
which complete data were received (solid line); com- 
pared with this is the approximation of an exponential 
time -dependence (dashed line). 

apparent rotation calculated using the methods 


to be described is shown in Figure 2. This plot and others like it for subsequent shadow passes 
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Figure 2— Change of apparent azimuth of magnetic field 
observed (solid line) during same shadow pass as in Fig- 
ure 1; compared is the predicted change (dashed line) 
calculated by assuming that the magnetic field did not 
change its direction. 


tered once every four sequences 
tion of the spacecraft about the spin axis at any time. 


attest to the accuracy of the calculations, 
using a time-dependent spin period of the form 
in Equation 1. 

To help clarify the method of analyzing the 
spacecraft aspect in the shadow, a brief de- 
scription of the pertinent facts of spacecraft 
operation is necessary. The direction of the 
sun relative to the spacecraft coordinate sys- 
tem is measured once every four telemetry 
sequences (1 sequence = 81.808 seconds). 

This is measured by counting the time from a 
fixed reference point in the telemetry sequence 
to the next "see sun" pulse from the onboard 
sun sensor.* This time (referred to as the 
"sun time") and the time between consecutive 
sun pulses (the raw spin period) are teleme- 
This enables one to calculate accurately the azimuthal orienta- 


Such an accurate determination of the spacecraft azimuth requires the calculation of a refined 
spin period. This is done by calculating the time between successive sun times and dividing by the 
integral number of revolutions made by the spacecraft. The number of revolutions N is calculated 
from the raw spin period. Thus the refined spin period is given by 


+ c - c 


( 2 ) 


where 

C 0 = the time between successive sun reference times, 

Cj = sun time at beginning of interval, 
c 2 = sun time at end of interval, and 

N = integral number of revolutions between sun sightings. 
N in turn is given by 


N ± A 



(3) 


where S is the raw spin period, and A is a fraction. The fraction A will be less than 0.5 as long 
as the time C 0 corresponds to less than approximately 40 sequences. Thus N is found easily by 
rounding the quantity (N + A) to the nearest integer. 


*AIMP (IMP-D) Technical Summary Description, GSFC Document X-724-67-87, March 1967. 
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If the spin period remained constant through the shadow, it would be a simple matter to deter- 
mine accurately the azimuthal orientation of the spacecraft. Actually the spacecraft spins up in 
the shadow; and, by the end of the shadow, a large error (of the order of 360 degrees) has accumu- 
lated in the satellite orientation computed from the optical aspect data. It is this error that must 
be corrected. 

The time when the spacecraft enters the shadow is measured to the nearest half sequence, 
except every fourth sequence when it is accurate only to the nearest sequence. This information 
is telemetered in a sun flag, which indicates whether the sun sensor is seeing the sun or not. 

Thus the situation is as follows: The spin period and sun direction are known accurately at 
the time of the last sun sighting before the shadow. Then there follow 0 to 4 sequences of meas- 
urements before the satellite enters the shadow 
and the spin period begins to change. Upon en- 
tering the shadow, the spacecraft begins to spin 
up, and the spin period begins to decrease. This 
decrease continues until the spacecraft again 
enters the sunlight whereupon it begins to de- 
spin again; and the spin period begins to in- 
crease, finally approaching its preshadow value. 

After leaving the shadow, there will be 0 to 4 
sequences before true solar aspect data are ob- 
tained. The total fractional change in the spin 
period is approximately 0.0015 (0.15 percent) 
over the period of the shadow— typically 40 min- 
utes (35 sequences). 

With this understanding of the problem, 
calculating the three free parameters in Equa- 
tion 1 is fairly straightforward. The following 
parameters are used (Figure 3 illustrates their significance): 

At. = time between the last true sun pulse and the beginning of the shadow, 

At = duration of the shadow, 

At f = time between shadow end and next true sun time, 
t = spin period before shadow, 

T = first spin period measured after shadow, 
t 2 = second spin period measured after shadow, 

C. = last sun time before shadow, 

C f = first sun time after shadow, and 
S = raw spin period. 



Figure 3— Spin period versus time for the same shadow 
pass as in Figure I;the solid part of the curve is derived 
from measured spin periods, and the dashed portion in the 
shadow is the approximation of Equation 1 with appro- 
priate values of the constants a, b, and t substituted. 
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In addition the following parameters will be used: 


r f = spin period at the end of the shadow, 

R(t) = actual number of revolutions of the spacecraft since last true sun before shadow 
(time "zero" on Figure 3), 

R 0 = R(t = At. + At s + At f ) = total number of spins made between true sun pulses, and 
N 0 = integral number of rotations of spacecraft between true sun pulses. 

Convenient units are as follows: 

At T s and r f s in seconds, 

C., c f , and S in 800-Hz clock counts, 

R, R 0 , and n in revolutions. 

It is often useful to express the At T s in units of sequence length (81.808 seconds); t 0 also will be 
measured in units of sequence length. 

COMPUTATION OF SATELLITE'S ROTATIONS 

The problem now becomes one of fitting a spin period of the form in Equation 1 to the satellite 
data. There are three free parameters in one; three independent equations are required to deter- 
mine these parameters. Two are obvious, stating that the spin period before and after the shadow 
must be equal to that measured. The third states that the integrated number of revolutions between 
the true sun times be equal to that achieved by the spacecraft. To this end, r f and R 0 must be 
determined since they are not measured directly. 

First, r f may be found by using a linear interpolation from r x and at the end of the 

shadow. (Earlier work also fits these data to an exponential to find r f ; but this was found to be 
unnecessary, the linear approximation being more than accurate enough for the short interval in- 
volved.) Thus (for At f in units of the sequence length) 


r 


f 



A little thought will reveal that R 0 is given by 


R 


o 



( 2 ) 


( 3 ) 


where N 0 is a constant to be determined; it may be found by integration if the spin period is known 
in terms of time. In this work, the spin period is assumed to be given by 


4 



t < At. 
1 


T. , 


r(t) = <a + be"^* ‘ At i )/t ol, At. < t < Lj , 


r f + ( T2 4 Tl ~) (t " At s “ At i) ' L 1 < 4 * L 2 


(4) 


where = At. + At s , and l 2 = L x + At f . Now R 0 may be obtained from: 


R 


o 



dt 

r(t) 


At i At S *0 . / 


a + be 


‘ At *Ar 


+ b 


At . 


( T 2 “ ^7) At f/ 4_ ' 


(5) 


Since initially the values of a, b, and t 0 are unknown, we still do not know the value of R 0 . 
However, the present goal is to obtain N 0 ; therefore, it is sufficient to know R 0 within ±0.5 revo- 
lution. The following approximate values of a, b, and t 0 give R 0 within ±0.1 revolution: 


a £ t f - 0.0001 second , 

b = r. - a, ^ (t. - T x j + 0.0001 second, and 
t Q = 24 sequences ^ 1963 seconds . 

Thus by rounding R 0 to the nearest integer, N 0 is determined. This is exact if the estimate 
of R 0 was within ±0.5 revolution; then Equation 3 may be used to determine R 0 accurately. The 
accuracy of N 0 is easily verified since an error shows up in the magnetic field data as an apparent 
rotation of one or more complete revolutions. When the foregoing information has been derived, 
it is a fairly straightforward matter to establish the three conditions that determine the three un- 
known parameters in Equation 1: a, b, and t 0 . 

The first two equations are obtained by requiring that the spin period be continuous at the 
beginning and end of the shadow. When the satellite enters the shadow, t = At. and r = r. so that 

a + b = t. . (6) 

At the end of the shadow, t = L x and r = r f . Thus 

a + be" ( At s^ t o) = T f . (7) 


The third and final condition is established by requiring that the predicted number of revo- 
lutions between the last real-sun pulse before the shadow and the first real-sun pulse after the 
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I III 


shadow is equal to the actual number of revolutions, r q . This requirement is expressed by Equa- 
tion 5. Since some of the terms in Equation 5 involve only the known parameters (i.e., do not 
involve a, b, or t 0 ), it is useful to define another quantity R s , the number of revolutions the space- 
craft makes during the actual shadow: 


R 


s 


At. 



At f 

T f ( T 2 " T l) At f/ 4 * 


Then the third condition may be rewritten as 


R 


s 




a + be ~ ( A * s / t o) 


a + b 


( 8 ) 


(9a) 


Since the change in the spin period during the shadow is observed to be only approximately 
1/1000 of the spin period, it is easily seen that b/a ~ 0.001 « 1. Thus it simplifies the problem 
to expand the logarithmic term in Equation 9a in powers of b/a. This results in 


R 


s 




1 - e ( At s //t o) , 


(9b) 


where only terms up to the first power of b/a have been retained. 

The next step is to solve Equations 6, 7, and 9b simultaneously for a, b, and t Q . From 
Equation 9b, 


~ e ~ ( At s/t °) 



( 10 ) 


From Equation 7, 


t o b 


( Ats/t °) - 77 ( l f ' a ) 


and from Equation 6, b = r. - a. Substituting the last two expressions into Equation 10 yields 


At s t Q 

R s -~r + tt ( r i ■ T f) = 0 • 


Upon rearrangement, this expression can be written 


a2 At 


- a + 


AT ( T i - r f) = 0 
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From Equations 6 and 7 it is seen that 


e“ ( At s / o) 



r . 


a 


a 


so that 


or 



Substituting the foregoing expression in Equation 10 yields 


R 



At 


- a + 



(ID 


( 12 ) 


It is possible to find an approximate solution to Equation 12 by use of Newton T s method.* This 
is an iterative method for finding the roots of an equation of the form f (X) = 0 . In essence, an 
initial estimate of the root X 0 is used to find an improved estimate x t from the equation 


Thus, in the present problem, 


1 df(X) 


Then 


f(a) 


0, where f(a) 


a 2 R 


At 


s 




When a has been determined this way, it is a simple matter to use Equation 6 to determine b, 
and Equation 11 to obtain t 0 . 


*Philip Franklin, "Methods of Advanced Calculus/* pp 26-28, New York: McGraw-Hill, 1944. 
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Integration of Equation 1 yields the actual number of rotations of the satellite; i.e., 


r 


t 


T . 


t < At . 


K(t) = <R, + 


t - At. t 0 b 


(t - At ,j 
1 - exp - - — j 

L n 


At £ < t < L, 


t - L, 


R 2 + 


r f + 


L ) (t - Lj 


L x < t < L 2 


where R 2 = At .yv. , and 


D _ D At s ^ 
R 2 R i + a ~ a 2 


At 

1 - exp j- — 


NUMERICAL EVALUATION OF PARAMETERS AND THEIR APPLICATION 

An application of the method previously described to the pass shown in Figure 2 is useful for 
illustrative purposes. The input parameters are as follows: 

r. = 2.2935 seconds, 
r f = 2.2909 seconds, 

At s = 30.5 sequences, and 
R = 1088.675 revolutions. 

s 

The resulting values for a, b, and t 0 are 

a = 2.2898 seconds, 
b = 0.003641 second, and 
t 0 = 23.96 sequences. 

With the additional information that 

r x = 2.29125, 
r 2 = 2.29151, 

At 4 = 1.5, 

At f = 4.0, 

and a knowledge of the actual spin period, the correction angles can be calculated as follows. If 


R u (t) 


J 0 
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where r u is the spin period used in the initial 
analysis, then the angle between the apparent 
sun direction (represented by the pseudo-sun 
pulse) and the true sun direction is given by 

t) = [R(t) " R u (t)] rotations 

= 360 [R( t ) ~ R u (t)] degrees. 

The correction angles for the pass shown 
in Figure 3 are listed in Table 1. If it is as- 
sumed that the magnetic field remained con- 
stant in direction during this shadow pass, then 
the apparent direction would have rotated 
through the angle i//( t). It is on this assumption 
that Figure 2 was plotted. The solid curve 
shows the azimuthal angle of the interplanetary 
magnetic field obtained using the pseudo- sun 
pulse for orientation. The dotted curve shows 
the relative change K t). 

In conclusion, the analysis presented here 
is an effective method for correcting the azi- 
muthal orientation of the Explorer 35 spacecraft 
during its passes behind the moon, where it can 
no longer use the sun as a reference. The ac- 
curacy of the method is better than ±10 degrees, 
but can only be checked empirically if it is as- 
sumed that the magnetic field does not change 
direction. The largest error would occur dur- 
ing the shadow since, at the beginning and at the 
end, the orientation has been fitted to the true 
direction obtained by real- sun sightings. 


Table 1 


Correction Angles for Pass Shown in Figure 3. 


t 

R(t) 

K 0 (t) 

'A(t) 

(sequences) 

(revolutions) 

(revolutions) 

(degrees) 

0 

0.0 

0.0 

0 

1 

35.6697 

35.6608 

3 

2 

71.3395 

71.3216 

6 

3 

107.0114 

106.9823 

10 

4 

142.6856 

142.6431 

15 

5 

178.3618 

178.2982 

23 

6 

214.0401 

213.9532 

31 

7 

249.7203 

249.6083 

40 

8 

285.4023 

285.2633 

50 

9 

321.0862 

320.9184 

60 

10 

356.7717 

356.5734 

71 

11 

392.4592 

392.2284 

83 

12 

428.1482 

427.8835 

95 

13 

463.8386 

463.5385 

108 

14 

499.5305 

499.1936 

121 

15 

535.2241 

534.8486 

135 

16 

570.9187 

570.5037 

149 

17 

606.6147 

606.1587 

164 

18 

642.3118 

641.8138 

179 

19 

678.0103 

677.4688 

195 

20 

713.7100 

713.1238 

211 

21 

749.4106 

748.7789 

227 

22 

785.1123 

784.4339 

244 

23 

820.8149 

820.0890 

216 

24 

856.5188 

855.7440 

279 

25 

892.2234 

891.3991 

297 

26 

927.9287 

927.0541 

315 

27 

963.6350 

962.7092 

333 

28 

999.3425 

998.3642 

352 

29 

1035.0503 

1034.0192 

371 

30 

1070.7590 

1069.6743 

391 

31 

1106.4687 

1105.3293 

410 

32 

1142.1787 

1140.9844 

430 

33 

1177.8889 

1176.7318 

417 

34 

1213.5979 

1212.4792 

403 

35 

1249.3059 

1248.2266 

389 
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